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Summary 

In order to achieve a better understanding and comparison of conventional and vertical 
junction laser-photovoltaic energy converters, we have constructed three-dimensional analytical 
models of their operation. These models simulate the conversion of laser power into electrical 
power and are constructed for both planar conventional and vertical junction photoconverters. 
The models are a set of linear partial differential equations and take into account many 
photoconverter design parameters. From the models, we extract the spectral response and 
current-voltage relations from solutions to the partial differential operators in terms of Green’s 
functions. The model is applied to Si photoconverters; however, the model may be used with 
other semiconductors. The conversion efficiency of the Si vertical junction photoconverter is 
47 percent, whereas the efficiency for the conventional Si planar device is only 17 percent. 
A parametric study of the vertical junction Si converter is then done in order to describe the 
optimum converter for use with a 1.06-//m Nd laser. The optimum Si vertical junction converter 
for a single junction is 1 mm deep by 1 cm wide by 0.035 mm thick. The junction is located 
2 /mi from one edge. The doping density for both acceptors and donors is 1 x 10 18 cm -3 . A 
reflectivity of 5 percent is assumed. The efficiency for this optimized vertical junction converter 
is 44 percent at 1 kW-cm -2 . 

Introduction 

Photovoltaic converters have many applications for satellite and space station power needs. 
In this paper, we are concerned with the application of photovoltaics for the conversion of 
laser power to electrical power in a space-based laser power system (ref. 1). We envision laser 
radiation on the order of 1.0 kW-cm -2 (ref. 1) at a specific wavelength that depends upon the 
type of laser considered. Semiconductors for use as laser power converters must be chosen such 
that their bandgap energy is less than the energy of the laser photons (ref. 2). In this study, we 
develop three-dimensional analytical models to compare planar conventional photoconverters 
with vertical junction photoconverters in order to achieve a comparison of how these energy 
converters operate. The photoconverter considered is Si which has a bandgap energy of 1.11 eV . 
For this type of photoconverter, we consider a Nd laser which emits 1.06 /tm radiation which 
is near the Si bandgap. The models for the conversion of laser power to electrical power are 
based upon a set of linear partial differential equations with appropriate boundary conditions. 
Solutions of these equations are obtained in terms of Green’s functions (ref. 3) for the regions 
n and p. From these solutions, various output quantities are calculated. In particular, we 
calculate the efficiencies of the converters, the spectral response, current voltage relationships, 
and contour plots of the current density at the junction edge. Temperature gradients are 
neglected. The method used is very general and will handle a wide variety of source terms. 

Model 

General 

The geometry for the three-dimensional models of an n/p photoconverter is illustrated in 
figure 1. We denote the region n by 

Rn = {x,y,z | |x| < a, |y| < 6, 0 < z < zj } 

and denote the region p by 

Rp = {x,y,z | |x| < a, |j/| < b, zj + u <z < H) 

The junction doping is assumed to be linear and abrupt at the junction edge. For the planar 
conventional photoconverter, we let 


g(z) = 0 O (1 - R)ot exp(-az) 


( 1 ) 



denote the generation rate of the hole-electron pairs at a depth z when light is incident 
perpendicular to the surface z = 0. We let g{a — x) denote the generation rate for the vertical 
junction photoconverters when light is incident upon the face x = a. In equation (1), <t>o in 
cm -2 -s -1 is the rate of incident photons per second per unit area, R is the fraction reflected 
from the surface, and a is the absorption coefficient in cm -1 at the wavelength considered. 

Conventional Photoconverter 

We consider first the model for a conventional planar photoconverter with light incident 
upon the face z = 0. For the region Rn , the diffusion of the excess carriers p = p n — p nQ in the 
n-material is governed by the equation 


H(p',Lp) = —g(z)/Dp ( x,y,z € Rn) 

where 

H(u; L p ) = v 2 u — £^2 ( 2 ) 

which is subject to the following boundary conditions: 


B^Pi x , Dp) BoP 
B(p; x; D p ) = -S 0 p 
B[p',y,D p ) = Sop 
B(p;y;D p ) = -S 0 p 
B(p; z; D p ) = -S p p 

p = 0 


(x = a) ' 
(x = —a) 

(y = b) 

( y = -b ) 

(z = 0) 

{Z = Zj) ; 


( 3 ) 


with 

u( \ du 

B(u;x; c) = -c— 

a boundary operator, and where c is a dummy parameter within the operator. The parameters 
in equations (2) and (3) are So, the surface recombination velocity at the sides; S p , the surface 
recombination velocity at the surface z = 0; D p , the diffusion coefficient of holes; and L p , the 
diffusion length of holes in an average lifetime. 

For the region i? p , the diffusion of the excess minority carriers n = n p — n PQ in the p-material 
is governed by the equation 


H(rv,L n ) = (x,y,z,ERp) (4) 

JJn 

which is subject to the following boundary conditions: 


B(n;x;D n ) = Sin 

(x = a) 

B(n;x-,D n ) = —Sin 

(x = —a) 

B(n;y;D n ) = Sin 

(y = b) 

B(n;y,D n ) = -Sin 

(y = -b) 

B(n; 2 ; D n ) = S n n 

( z = H) 

n = 0 

(z = Zj + w) 


with parameters Si, the surface recombination velocity at the sides; S n , the surface recombi- 
nation velocity at the surface z = H\ D n , the diffusion coefficient of electrons; and L n , the 
diffusion length of electrons. 

Associated with the boundary value problem using equations (2) and (3) is a Green’s 
function G\(x, y, z; xq, yo, zq) which can be obtained by the method of separation of variables. 
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This enables us to represent the solution to the boundary value problem of equations (2) and (3) 
in the form: 


P — Pn ~ Pno — 



-g(zp) 

Dp 


G \ ( x , y, z\ XQ,yo, zq ) dx o dyo dzo 


(6) 


Similarly, there is a Green’s function G2{x,y,z;xo,yo,zo) associated with the boundary value 
problem using equations (4) and (5) which enables us to represent the solution to equations (4) 
and (5) in the form: 


rH rb ra —g( z q) 

n = n p - n P0 = / / / G 2 (x, y, z\ x 0 , y 0 , zq) dx Q dy 0 dz 0 

Jzj+u J—b J—a 


(7) 


The functions G i and G2 can be expressed as a double Fourier series in terms of the 
eigenfunctions associated with the Sturm-Liouville problem (ref. 4): 


d 2 u 
1 r 2 + 


A 2 u = 0 


(8) 


with the following boundary conditions: 

B{u\ r; D ) = S311 (r = h) 
B(u ; r; D ) = S3U (r = — h ) 

The eigenfunctions for this problem are given by 


u n = «(r; D, S3, h, A n ) = cos A n r 


(9) 


where the eigenvalues A n satisfy the equation 

tan = — — ($h = 

Note that for large values of n, f n+ i approaches mr. Also the eigenfunctions are orthogonal 
over the interval (—A, h) with norm squared given by 

Kll 2 = h+ COs2 2 (<Tn = A n h) (10) 

We can then define the eigenfimction associates with the regions R n and R p as 

U m {x) = u(x; D p , Sq, a, X m ) = cos A m x 


V m {y) = u(y, D p , S 0 , b, 77m) = cos r) m y 
U^x) = U(x;D n ,Si,a,\n) = cos X* m x 


( 11 ) 


lm(g) = Si, 6, r/m) = cos r]* m y 


The Green’s functions G\ and G 2 can be represented in the form of a double Fourier series 
symmetric in x and y given by 


00 00 

<?i = E£ 


1=1 i=l 


'hjj Ujjxo) Uj{x) Vj{yo) Vj{y) 

A i WUjlfm 2 


( 12 ) 
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where 


i=lj=l 2 


^..-IY 1 (z 0 )Y 2 (z) 

* U 3 -\Y 1 {.z)Y 2 {zq) 


thijUjMUjWVfiv o)v?(v) 


\ u n m 


(0 < ZQ < z) 
(z < ZQ < Zj) 


l/)n ■ ■= I r 3(*o) **(*) ^ + “< z 0< z ) 

^ \ ^3(4 Vito) {z<ZQ< H) 


7 ij sinh {zj/'iij) 


Y 4 (zi + u) 


Yi(z) = cosh ( — ) + sinh (— 


14W-^ nh 

sinh {z/'iij) 


Ys(z) = sinh 


z — ( Zj + u) 


Y 4 (z) = cosh (— — - | sinh (*L — £ 


having the parameters 7 y and ay defined by 


7?.= 

nj 


*ij = 


1 + L^{ A? + r/|) 
i + LKxf + nf) 


From equations (6) and (7), we can obtain the current densities as follows: 


Jp — —qDn 


d ( Pn ~ Pnp ) 


(z = Zj) 


t _ _r> d ( n P~ n Po) , _ , i 

Jn — qD n — ( z — Zj + u>) 


This gives the total short circuit current J sc = J p + J n - 
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By an analysis of Green’s second formula, we can calculate the dark currents associated 
with the regions R n and Rp. In the region Rn, we obtain a modification of the boundary value 
problem using equations (2) and (3). The modification being that in equation (2), we let g = 0, 
and in equation (3), we require that 


Pn-Pn 0 =Pn 0 exp 


(SH 


{z = Zj) 


where Vj is the junction voltage. This boundary value problem has the solution 


, *- p Wl£M exp (w)- 


dG i 


dz o 


dxQ rfj/o 


Z 0 =Zj 


from which we can obtain the dark current opposing the photocurrent as 

Jpd = ~qDp grad (p n - p„ 0 ) • fc {z = Zj ) 


(17) 


( 18 ) 


Similarly, in the region Rp, we obtain for the model of the dark current a modification of 
the boundary value problem with equations (4) and (5). The modification being that in 
equation (4), g = 0, and in equations (5), we require that 


(§)- 


n p “ n P0 = [ ex P 
This boundary value problem has the solution 

rb 


(Z ~ Zj + w) 


Tin — n 


PO 


r r (qVj\ 9g 2 

= L L nw h \kt) ~ 1 


ZQ=Zj+UJ 


dx o dyQ 


which, in turn, produces the dark current 

Jnd = ~qD n grad (n p - n Po ) ■ (- k ) (z = Zj + w) 


(19) 


( 20 ) 


The total dark current is then a summation of the results in equations (18) and (20) and we 
obtain 


d d, — Jpd 4" Jnd ~ Jo ex P ^ ^ 


( 21 ) 


We use the recombination current density J Tec from reference 4 and calculate the total current 
by integration with respect to x and y over a plane represented by z = Constant. 

The current-voltage relation is obtained from the equation relating the current and voltage 
for the equivalent circuit diagram of figure 2 and we find 


I ( 1+ R sh ) hc R 9h Id 


where R s is the series resistance and R s h is the shunt resistance. 

Vertical Junction Photoconverter 

For the vertical junction photoconverter, we assume that light is incident upon the face 
x = a. In the boundary value problems using equations (2), (3), (4), and (5), the generation 
term g(z ) is replaced by g(a — x). Also in equations (6) and (7), the form of the double Fourier 
series representing the Green’s function must be changed because in the vertical junction case 
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we no longer have symmetry in the x-direction. We consider the following Sturm-Liouville 
problems (SLP) in the regions Rn and Rp. In the region Rn, we consider 


= W" + r] 2 W = 0 (0 < z < Zj) 

subject to the boundary conditions 

W = 0 (z = Zj ) 

B(W;z,D p ) = —S P W (z = 0) 

In the region Rp we consider the SLP 


( 22 ) 


L(W*) = W*" + r, 2 W* = 0 (z j + u<z<h) 
subject to the boundary conditions 

W* = 0 (z = zj + uj) 
B{W*-,z,D n ) = S n W* ( z = H ) 

The eigenvalues of the SLP in equation (22) satisfy 

-D p (r)iZj) 

— rr- = tan ( t = 1, 2, 3, . . .) 

Z)pZj 

and the eigenvalues of the SLP in equation (23) satisfy 

-D n Mti) ± , , 

— ^ — = tan nfl ( i = 1, 2, 3, . . .) 


(23) 


(24) 


(25) 


with b! — H — Zj - uj. The corresponding eigenfunctions for the regions Rn and Rp are, 
respectively, 


Wi(z) = sin r)i( Z j - z ) (i = 1, 2, 3, . . .) 


and 


w i (*) = sin r)i (z - zj - u) (i = 1, 2, 3, . . .) 
which have the square norms 


1 sin 2 rjiZj 
- 2 Z 1 ~ 4 m 


(*' = 1, 2, 3, ...) 


, ra2= l^_lsin|^ ( j = i, 2, 3, . . .) 

•I 


(26) 

(27) 

(28) 
(29) 


For the vertical junction case, the Green’s functions for the regions Rn and Rp are represented 
as the double Fourier series 




oo oo 

92(x,y,z;x 0 ,yo,z 0 ) = EE 

i=l j=l 


Vito) wjto) v i*(vo) wjizo) ifivj 
Iiv-i 2 i|w7ii%. 


(30) 

(31) 
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where Vf and V{ have been previously defined in equation (11). The remaining terms in 
equations (30) and (31) are determined from the relations: 


\ -Yi(x) y 2 (x 0 ) 

(-a <xq<x) 

ij = < 


( -yi(zo) Y 2 (x) 

(x < xo < a) 

( -y 3 (x) y 4 (x 0 ) 

(—a < xo < x) 

Tp4ij = { 


{ -y 3 (x o) y a (x) 

(x < xq < a) 


Yi{x) = Fi[x\Tij,So,D p ) 

Y 2 (x) = n(-z) 
y 3 (x) = Fi(x; £ij, Si,D n ) 

Y 4 (x) = Y 3 (-x) 

Fi(x;A,B,C) = cosh (^p) + ~§~ sinh (~p) 


1-2 

2 

13 1 + Lp + tf) 

2 L\ 

£ij ~ 1 + Ll^f + rjf) 
Pij — ^2 (Fiji , S0j D p ) 

Pij ~ ( e ij 5 SliD n ) 


( 1 B 2 A\ . ^ (2a \ 2 B ^ (2a\ 

{a + -&) smh (x) + -c C 0 Sh {A> 


f 2 (a,b,o)=i j + - 51 

In equations (2) and (3), we replace g(z) by g(a — x), then in the region Rn the density of 
excess minority carriers can be represented in the form 

fZj rb ra q ( q , — X{\ 1 

Pn Pn o = / / / 91 (®j Vi VOi ^o) dxo dy q dzQ 

Jo J-b J-a Bp 

Similarly in equations (4) and (5), we replace g(z) by g(a — x ) and find that the solution for 
the excess minority carriers in the region Rp can be represented in the form 


Tip TlpQ — 



g(q - xp) 


D 


n 


g2(x,y)Z]yo,zo) dxo dy 0 dzQ 


With the photocurrent relations of equation (18), we can produce the photocurrent responses 
at the junction edges Zj and zj + w from which we can calculate the total short circuit current. 

Computer programs simulating both the conventional and vertical junction photoconverter 
responses have been developed. As a check on the results obtained from these three-dimensional 
models, we constructed one-dimensional models of both the conventional and vertical junction 
photoconverters, and the results from that study can be found in reference 5. 
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Results and Discussion 


Comparison of Conventional and Vertical Junction Converters 

Figure 3 gives the efficiencies of a conventional converter and a vertical junction converter for 
Nd laser radiation. These efficiencies were calculated by using the nominal parameters shown 
in table I. The efficiency of the vertical junction converter at an input power of 1 kW-cm” 2 is 
47 percent, whereas under the same conditions, the efficiency of the conventional converter is 
only 17 percent. Because of the higher efficiency of the vertical junction converter, a detailed 
parametric study has been done for this converter. 

Parametric Study of Vertical Junction Converters 

Since the efficiency of the vertical junction converter is significantly higher than the efficiency 
of the conventional converter, a parametric study of the vertical junction converter has been 
done. Baseline parameters shown in table II were assigned individually to determine the 
optimum, achievable value for each parameter. The best set of these values was then assigned 
to describe the optimum vertical junction converter for use with a 1.06-//m Nd laser. 

Application of Model to Si Vertical Junction Converter Nd Laser System 

It has been shown previously that Si vertical junction converters are appropriate for use 
with a Nd laser to form a space power system (refs. 2 and 5). We have applied our three- 
dimensional model for vertical junction converters to such a system. The effects of various 
design parameters are shown, and the optimum single junction converter is described. 

The set of converter baseline parameters chosen is shown in table II. These parameters were 
varied individually in order to determine their optimum values and then assembled to describe 
the optimum single junction converter. Each parameter is discussed separately. 

Carrier Concentration 

Figure 4 shows the effect of carrier concentration on converter efficiency for our baseline 
conditions. The peak baseline efficiency of 47.5 percent occurs at an acceptor concentration 
of 1 x 10 18 carriers-cm” 3 . The peak baseline efficiency of 51.2 percent occurs at a donor 
concentration of 1 x 10 18 carriers-cm -3 . We chose for our optimum carrier concentrations both 
donor and acceptor concentrations of 1 x 10 18 carriers-cm” 3 . 

Converter depth . Figure 5 shows the effect of converter depth on converter efficiency. The 
efficiency increases from 20.5 percent at a depth of 5 pm to 47.8 percent at a depth of 35 /zm. 
We have taken 35 pm as the optimum converter depth. 

Converter thickness . Figure 6 shows the effect of converter thickness on converter efficiency. 
The efficiency increases to a maximum baseline efficiency of 45.5 percent at a converter thickness 
of 1000 pm . This thickness of 1000 pm is used as the thickness of our optimum converter. 

Junction position • Figure 7 shows the effect of junction position on converter efficiency. As 
the depth of the n-layer (junction position) decreases, the converter efficiency increases. From 
a fabrication consideration, we have chosen 2 pm as a practical, achievable junction position 
for our optimum converter. The peak baseline efficiency is 51.2 percent. 

Series resistance . Figure 8 shows the effect of series resistance on the converter efficiency. 
The efficiency decreases from 45.5 percent at 0.1 mO to 7.7 percent at a series resistance of 
1 ohm. For our optimum converter, we have chosen a series resistance of 1 mfi corresponding 
to 20 pm of silicon thickness. The baseline efficiency with this series resistance is 45.4 percent. 

Temperature . Figure 9 shows the converter efficiency as a function of temperature. Note 
that the efficiency decreases with increasing temperature. The efficiency at our operating 
temperature of 300 K is 45.4 percent. Operation of Si above 500 K may not be practical. 
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Surface recombination velocity. Figure 10 shows the effect of n-surface recombination velocity 
and p-surface recombination velocity on converter efficiency. As the n-surface recombination 
velocity increases from 0.01 m-sec -1 to 100 m-sec -1 , the efficiency decreases from 49.4 per- 
cent to 31.1 percent. As the p-surface recombination velocity increases from 0.01 m-sec to 
100 m-sec -1 , the efficiency decreases from 47.6 percent to 38.7 percent. According to refer- 
ence 6, a surface recombination velocity of 1 m-sec -1 is achievable in Si; therefore, we have 
chosen 1 m-sec — 1 as the recombination velocity of our optimum, achievable converter. The 
baseline efficiency at this surface recombination velocity is 45.5 percent. 

Optimum converter. From this parametric study, a set of optimum parameters has 
been determined. These are shown in table III. The efficiency of the optimum converter is 
43.9 percent at 1 kW-cm -2 . Specific results from this computer study were presented in 
graphical form to illustrate the importance of the parameter changes. Major factors to consider 
in the design of an efficient vertical junction, high-intensity photoconverter are the surface 
recombination velocity and the series resistance. Figure 11 shows the converter efficiency as 
a function of input power for the optimized vertical junction converter. The efficiency varies 
from 33 percent at 1 kW-cm -2 to 44 percent at 1 kW-cm -2 . 

Concluding Remarks 

We have presented a three-dimensional model for the analysis of a conventional and a 
vertical junction photoconverter for the conversion of laser power into electrical power. The 
model requires the representation of Green’s functions in the form of a double Fourier series. A 
set of nominal values was assigned to the parameters in the model, and then a comparison 
of the conventional and vertical junction photoconverter was undertaken. The calculated 
efficiency of the vertical junction converter is 47 percent, whereas the calculated efficiency 
of the conventional converter is 17 percent. A parametric study of the higher efficiency vertical 
junction converter has been done. The efficiency at 1 kW-cm 2 for our optimized converter is 
44 percent. 

NASA Langley Research Center 
Hampton, Virginia 23665-5225 
May 3, 1988 
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Table I. Nominal Values Used in Initial Calculations 


Parameter 

Conventional 

Vertical junction 

Dimensions (see fig. 1): 



Thickness, a, mm 

5 

5 

Width, 6, mm 

Junction depth, zj, fim 

5 

5 

5 

20 

Mobility for — 


Region n, cm 2 -V -1 -s -1 

711.3 

711.3 

Region p, cm 2 -V -1 -s -1 

89.7 

89.7 

Laser wavelength, // m 

1.06 

1.06 

Surface recombination velocity on — 

Surface n, cm-s -1 

1.0 

1.0 

Surface p, cm-s" 1 

1.0 

1.0 

Series resistance, fi 

0 

o 

Shunt resistance, fi 

1.0 x 10 6 

1.0 x 10 6 

Doping density for — 


Na (acceptor), cm -3 

1.25 x 10 17 

1.25 x 10 17 

Nd (donor), cm" 3 

5.0 x 10 19 

5.0 x 10 19 

Lifetime for — 


Region n, s 

10.0 x 10" 4 

10.0 x 10 -4 

Region p, s 

5.0 x 10 -5 

5.0 x 10 -5 

Laser power, kW-cm~ 2 

1.0 

1.0 

Temperature, K 

293 

293 




Table II. Baseline Parameters Used for Vertical Junction Case Parametric Study 


Laser power, kW-cm -2 1 

Laser wavelength, /zm 1.06 

Dimensions (see fig. 1): 

Thickness, a, mm 0.65 

Width, b, mm 5 

Depth, H, mm 0.02 

Junction position, /um 10 

Surface recombination velocity on — 

Surface n, cm-s -1 1.0 x 10 3 

Surface p, cm-s -1 1.0 x 10 3 

Series resistance, fi 0 

Shunt resistance, fl 1.0 x 10 6 

Doping density for — 

Na (acceptor), cm -3 1.25 x 10 17 

Nd (donor), cm -3 5.0 x 10 19 

Mobility for — 

Region n, cm 2 -V -1 -s -1 711.3 

Region p, cm 2 -V -1 -s -1 89.7 

Lifetime for — 

Region n, s 1.0 x 10 -3 

Region p, s 5 x 10 -5 

Temperature, K 300 


Table III. Parameters for Optimum Converter 


Laser power, kW-cm -2 1 

Laser wavelength, j/m 1.06 

Dimensions (see fig. 1): 

Thickness, a, mm 1000 

Width, 6, mm 5 

Depth, H, mm 35 

Junction position, //m 20 

Surface recombination velocity on — 

Surface n, m-s -1 1 

Surface p, m-s -1 1 

Series resistance, Q 0.001 

Shunt resistance, fi lx 10 6 

Doping density for — 

Acceptors, cm -3 lx 10 18 

Donors, cm -3 lx 10 18 

Temperature, K 300 



Laser 




H 


y z 


(a) Conventional converter. 



Figure 1. Geometry for three-dimensional model. 


Efficiency, % 


R s 



Figure 2. Equivalent circuit diagram. 
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Figure 3. Efficiency versus input power for conventional and vertical junction converters. 
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Figure 4. Efficiency versus acceptor and donor concentration for vertical junction converter under 
baseline conditions. 


Baseline parameters 



Converter depth, pm 

Figure 5. Efficiency versus converter depth for vertical junction converter under baseline conditions. 
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Figure 6. Efficiency versus converter thickness under baseline conditions. 



Figure 7. Efficiency versus junction position for vertical junction converter under baseline conditions. 
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Series resistance, m £2 


Figure 8. Efficiency versus series resistance for vertical junction converter under baseline conditions. 



Temperature, K 


Figure 9. Efficiency versus temperature for a vertical junction converter under baseline conditions. 
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Figure 10. Efficiency versus n-surface recombination velocity and p-surface recombination velocity 
under baseline conditions. 
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